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The average unsigned magnetic flux density in magnetograms of the quiet Sun is generally dominated by instrumental 
noise. Due to the entirely different scaling behavior of the noise and the solar magnetic pattern it has been possible to 
determine the standard deviation of the Gaussian noise distribution and remove the noise contribution from the average 
unsigned flux density for the whole 15-yr SOHO/MDI data set and for a selection of SDO/HMI magnetograms. There 
is a very close correlation between the MDI disk-averaged unsigned vertical flux density and the sunspot number, 
and regression analysis gives a residual level of 2.7 G when the sunspot number is zero. The selected set of HMI 
magnetograms, which spans the most quiet phase of solar activity, has a lower limit of 3.0 G to the noise-corrected 
average flux density. These apparently cycle- independent levels may be identified as a basal flux density, which represents 
an upper limit to the possible flux contribution from a local dynamo, but not evidence for its existence. The 3.0 G HMI 
level, when scaled to the Hinode spatial resolution, translates to 3.5 G, which means that the much higher average flux 
densities always found by Hinode in quiet regions do not originate from a local dynamo. The contributions to the average 
unsigned flux density come almost exclusively from the extended wings of the probability density function (PDF) , also 
in the case of HMI magnetograms with only basal-level magnetic flux. These wings represent intermittent magnetic 
flux. As the global dynamo continually feeds flux into the small scales at a fast rate through turbulent shredding, a 
hypothetical local dynamo may only be relevant to the Sun if its rate of flux build-up can be competitive. While the 
global dynamo appears to dominate the magnetic energy spectrum at all the resolved spatial scales, there are indications 
from the observed Hanle depolarization in atomic lines that the local dynamo may dominate the spectrum at scales of 
order 1-10 km and below. 
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1. Introduction 

Dynamo processes are responsible for the gener ation of 
macr o scopic magnetic fields from feeble seed fields (jLarmod 
Il919t lElsassed Il946l Il956[) . The solar dynamo, which is 
the source of the 11-year activity cycle with its sunspots, 
fiares, CMEs, prominences, etc., is a prototype of a cos- 
mic dynamo. The basic ingredients of the Sun's dynamo 
are the same as for the galactic and planetary dynamos: 
interactions between magnetic fields and turbulence in an 
electrically conducting and rotating medium. The rotation 
breaks the left-right symmetry of the turbulence through 
the Coriolis force, thereby generating a large-scale net he- 
licity that is the s ource of the dynamo-produ ced magnetic 
field (iParkerilTQSa iSteenbeck fc Krausdl 19691) . 

While the symmetry breaking mechanism can produce 
magnetic patterns on large scales, turbulent shredding en- 
sures that the magnetic scale spectrum extends over many 
orders of magnitudes, down to the magnetic diffusion limit 
near 10-100 m, far below the observationally resolved scales 
on the Sun (cf. IStenfioll2012[ ). The magnetic energy that is 
injected at large scales by the dynamo quickly cascades all 
the way down the scale spectrum. The turbulent scale re- 
distribution of the globally generated magnetic energy gives 
the Sun's magnetic field a fractal-like appearance. 

In the context of numerical simulations of magneto- 
convection the concept of a small-scale, "local dynamo" 
that operates near the solar surface has often been used 
to account for much of the small-scale magnetic structuring 



(lKazantsevlll968l:lPetrovav fc SzakaIvl ll993l:[Cattaneolll999l: 
iBrandenburg fc Subramanianll2005ir in contrast to the pre- 
viously discussed "global dynamo" , the local dynamo does 
not need a rotating medium or the breaking of the left-right 
symmetry of the convective motions but simply depends on 
the way convection tangles the magnetic field. However, in 
the great majority of the numerical simulations done to 
date, the amount of magnetic energy that is generated at 
small scales depends on the initial conditions of the simula- 
tion, and therefore they do not represent any local dynamo. 

A local dynamo only exists in the simulations if, with- 
out input from any global dynamo, magnetic energy can 
be built up from an initial seed field, and if the results of 
the simulations are independent of the choice of seed field. 
One clear demonstration t hat such a local dynamo i s pos- 
sible on the Sun is that of IVoeler k Schiisslerl (|2007[ ). but, 
as the authors point out, their results depend on chosen 
simulation idealizations, like the subgrid model for the vis- 
cosity. Therefore this demonstration cannot be interpreted 
as proving that a local dynamo contributes significantly to 
the observed small-scale magnetic fields on the real Sun. 

The question how relevant a local dynamo is for the Sun 
cannot be determined by numerical simulations alone, be- 
cause their results depend on the idealizations and simplifi- 
cations that have to be introduced. Instead we have to look 
for empirical answers derived from the observed properties 
of the Sun, trying to separate the relative contributions 
from the global and local dynamos. 
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We know that the cyclic behavior of the Sun's magnetic 
field (22-yr magnetic cycle with the 11-yr sunspot cycle) 
is the result of global dynamo action. In contrast the local 
dynamo is related to the quiet Sun and produces magnetic 
fields that are statistically time invariant. Therefore a local 
dynamo can only contribute to a possible "basal flux", a 
minimum magnetic flux level that is always present every- 
where and every time on the Sun, even in the complete ab- 
sence of sunspots and solar activity. The observe d presence 
of a b asal X-ray flux for slow-rotating cool stars (jSchriived 
|1987( ) has been interpreted as empirical evidence for the ex- 
istence of dynamo action in non-rotating plasmas without 
symm etry breaking from the Coriolis effect (Bercik et al.l 

In the next sections we will use data fr om SOHO/MDI 
(IScherrer et al.l[r9 95) and SDO/HMI (Sch errer et al.ll20Tl 
ISchou et al. 2012 ) to determine the Sun's basal magnetic 
flux density level. Since the apparent, time-invariant flux 
level has its main contribution from instrumental noise, 
we have developed a technique to statistically separate the 
noise and solar contributions from each other, to find the 
noise-free, intrinsic basal flux level. Such a separation is 
possible because the noise and the solar field obey com- 
pletely different scaling laws. The basal flux density level is 
resolution dependent, but in a way that is well defined from 
a known scaling law. The low value of the determined basal 
flux density allows us to place a tight limit to the possible 
contribution of a local dynamo to the observed small-scale 
magnetic fields on the Sun. 



2. Correlation between unsigned flux and sunspot 
number 

Recently a global analysis of the complete set of 96 minute 
cadence SOHO/MDI full disk magnetograms, 73,838 of 
them covering the 15-yr period May 1996 - April 2011, 
was carried out with the a hn of exploring the properti es of 
bipolar magnetic regions (jStenflo fc KosovichevI [20121 ) . As 
a byproduct of this analysis we determined for each magne- 
togram the average value -Bavc of the unsigned vertical flux 
density \By\ over a circular region around disk center, for 
which the normalized radius vector r/r© was less than 0.1, 
0.2, . . ., 0.8, and 0.9, respectively, including all pixels in the 
average (regardless of whether they occur in sunspots or 
not). This allows us to explore how well the value of i?avc 
correlates with the sunspot number for the various sizes of 
the averaging regions. 

The vertical flux density By was obtained from the 
line-of-sight component _B|| of the magnetograms through 
the simple relation By = where is the cosine of 

the heliocentric angle of the respective pixel in the im- 
age. This procedure is based on the assumption that the 
fields are on average oriented in the vertical direction. While 
this is the only practically feasible assumption that can be 
made, it is physically justified (at least for statistical pur- 
poses outside active regions) by the circumstance that most 
(> 90%) of the net flux recorded with the MDI 4arcsec 
resolution comes from collapsed, kG type flux elements 
(|Howard fc Stenflol[T97l IStenflol Il973l [l99l . Such strong 
fields get pushed in the vertical direction by the buoyancy 
forces in the photospheric layers where the fields are mea- 
sured. Higher (above the observed layers) the field starts 
to become force-free and may increasingly deviate from the 
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Fig. 1. Time series of the average of the unsigned verti- 
cal fiux density, smoothed with a 1-month wide time win- 
dow (solid lines) . Left panel: averaging over the disk region 
r/vQ < 0.1. Right panel: averaging over the disk region 
r/vQ < 0.9. The dashed curves use a second-order polyno- 
mial in the sunspot number (also smoothed with a 1-month 
time window) to fit the solid curves. The fit gets tighter as 
the disk averaging area increases. 

vertical. Errors due to this assumption will contribute to 
statistical scatter. The circumstance that the relation be- 
tween the disk averaged flux density and the sunspot num- 
ber turns out to be so tight, as we will see next, lends 
support to this procedure. 

Figure[T]shows as the solid curves the value of Save, time 
smoothed with a 1-month wide window, for two choices 
of the averaging region: r/r© < 0.1 (left panel), and 
r/rQ < 0.9 (right panel). The dashed curves represent the 
monthly averages of the international sunspot number Rz , 
converted to a magnetic field scale B^^ through the second- 
order relation 



Br^ =bo + hRz + b2Ri 



(1) 



This relation is fitted to the solid curves with the coeffi- 
cients as the three free fit parameters, which gives us the 
dashed curves. 

We see from Fig. [T]that the disk-averaged value of \By\ 
(right panel) can be nearly perfectly modeled, the fit is 
surprisingly tight, implying that there is an almost one-to- 
one relation between i?avc and the sunspot number. If we 
however only use the innermost portion of the solar disk 
for the averaging of \By\ (left panel), then Save fluctuates 
rather wildly, in particular during times of high solar activ- 
ity, and the fit with the sunspot number is poor. The reason 
is of course that the magnetic field in the disk center re- 
gion is not representative of the overall level of magnetic 
activity on the Sun. The presence or absence of individual 
low-latitude active regions govern the behavior of near 
disk center. One needs to account for the cumulative mag- 
netic contributions over a large portion of the disk to get a 
good representation of solar activity. 

It is nevertheless surprising that the sunspot number R^ 
turns out to be such a nearly perfect index for the overall 
magnetic activity of the Sun as represented by the disk 
average of the unsigned vertical fiux density. The largest 
deviations between the model fit and the data in the right 
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panel of Fig. [T] occur in regions of significant data gaps 
in the MDI time series (1998.48 - 1998.89 and 1998.97 - 
1999.17), so the physically relevant fit is even better than 
it may appear in the plot. As the size of the averaging 
region is increased, the fit between the Rz model and i?avc 
becomes quite tight already for r/r© < 0.5, and the fit 
improves only slowly as we go to still larger averaging areas. 
For r/vQ < 0.9 (right panel of Fig. [IJ the fit parameters of 
the model have the values (in G) bo — 19.0, bi — 0.24, and 
62 = -5.5 X 10-4. 

Coefficient bo, as defined by Eq. ([T|), can be interpreted 
as representing the apparent basal fiux density, the min- 
imum level reached in the absence of sunspots. However, 
instrumental noise contributes to a fictitious non-zero basal 
level, since the noise does not average out to zero when we 
average over the absolute, unsigned flux densities. Recently 
a careful noise analysis of the MDI magnetograms by 
iLiu et ahl (|2012f ) showed that the instrumental noise is in- 
deed of the same order as our apparent basal flux density 
level of 19 G. This indicates that the intrinsic level of the 
basal flux density must be much smaller. A determination 
of its actual value requires a very careful separation and 
removal of the noise contribution. The procedure for doing 
this is described in the next section. 

3. Statistical removal of MDI noise 

Fortunately it is possible to statistically separate the noise 
contribution from the intrinsic solar contribution, since the 
Gaussian noise and the fractal-like solar magnetic pat- 
tern obey entirely different scaling laws. The two contri- 
butions however combine in a way that is neither linear 
nor quadratic. Before we can model the scaling behavior 
we need to determine the form of the non-linear relation 
between the noise-affected and noise-free average unsigned 
flux densities, which will be done in the next subsection. 

3.1. Convolution of noise with the intrinsic PDF 

Let P{By) be the area normalized PDF (probability density 
function) for the intrinsic, noise- free flux densities. Then the 
average unsigned flux density is 

Savc = J PiB,)\By\dB,. (2) 

As shown in a detailed anal ysis of Hinod e SOT/SP data 
for the quiet Sun disk center (jStenfldl2010[ ). P{By) is char- 
acterized by an extremely narrow core peak centered at 
By = 0, which can be modeled by a stretched exponential. 
This peak is surrounded by quadratically declining damp- 
ing wings that extend out to the kG region. The contribu- 
tion to Save is completely dominated by the contribution 
from the damping wings, since the inner PDF core is so 
narrow. 

The apparent, noise-affected PDF, Papp, is obtained 
from the intrinsic P through direct convolution with the 
noise distribution, which can safely be assumed to have 
a Gaussian shape. For a Gaussian distribution with stan- 
dard deviation cr, the average of the unsigned fleld value 
of the distribution is 0.798 cr. If the intrinsic PDF were 
also Gaussian, then the convolved distribution would re- 
main Gaussian with a sigma that is obtained by adding the 
two individual sigmas quadratically. In contrast, when con- 
volving two Lorentzian profiles (which have quadratically 
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Fig. 2. Relation between the noise-free average unsigned 
vertical flux density (vertical axis, in the text denoted Save) 
and the corresponding noise-affected flux density (horizon- 
tal axis, in the text denoted i?app), for two values of a for 
the noise distribution, 18 and 34 G. The solid curves are 
obtained through numerical convolution of PDFs that are 
representative of quiet-sun magnetic flelds (as derived from 
Hinode analysis) with Gaussian noise distributions, while 
the dashed curves are analytical representations of the solid 
curves in terms of Eqs. ([3]) and 

declining wings, like the solar PDFs), their half widths add 
linearly. In the present case, however, we are convolving a 
Gaussian with a function that is more Lorentzian like. We 
would therefore expect that the average unsigned fleld val- 
ues would combine in a way that is neither quadratic nor 
linear, but somewhere in between. In the following we will 
determine the form of this relation. 

The noise- affected average unsigned flux density i?app, 
which is obtained from Eq. ^ if we replace P by Papp, 
can be determined directly from any given value for the 
noise a by direct numerical Gaussian convolution of P. To 
obtain this relation also as a function of Pave, representing 
solar regions with varying amount of magnetic flux, we have 
used the analytical shape of the symmetrized intrin sic PDF 
that represented our Hinode disk center data set ()Stenflol 
i20100 . and produced a family of P(Pu) distributions by 
varying the value of the damping parameter that governs 
the strength of the damping wings of the PDF and there- 
fore also the value of Pave- The result of these numerical 
convolutions and integrations are shown by the solid curves 
in Fig.O computed for two chosen values of the noise cr, 18 
and 34 G. Note that the average unsigned field values for 
these two Gaussians are smaller than these sigmas by the 
factor 0.798 and are thus 14.4 and 27.1 G, respectively. 

Next we have replaced the Hinode-type shape of P 
with a Lorentzian function and repeated the calculations. 
Convolution of a Lorentzian function with a Gaussian gives 
the well-known Voigt function. The resulting curves are 
practically indistinguishable from the solid curves in Fig. [5] 
and are therefore not plotted. This result is not unexpected, 
since the non-Lorentzian core of the Hinode PDF P does 
not contribute significantly to the value of Pave, and the 
PDF wings have the same shape as a Lorentzian function. 
Therefore our modeling is insensitive to the true shape of 
the PDF core region. 
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As we will later use inversion techniques to determine 
the actual noise in the magnetograms and calculate the 
scaling law for the noise-affected flux density, we need to 
represent the numerical results in Fig. [5] (the solid curves) 
in an analytical form, valid for arbitrary choices of the noise 
tr. Based on the considerations that we have made before 
about the way in which Gaussian and Lorentzian functions 
combine, our analytical model for the conversion of the 
noise- affected B^pp to the noise-free Save for a given value 
of the standard deviation a of the noise distribution is 

i?avo=[Ba"pp - (0.798 an (3) 

For two Gaussian functions a would be 2.0, for two 
Lorentzian functions it would be 1.0. By choosing a 
(through trial and error) to be 

a = 1.36- 0.004 cr + 0.0034 Bapp, (4) 

we obtain the dashed curves in Fig. [51 which can be seen 
to be nearly perfect representations of the exact, numerical 
relations. More importantly, our analytical model remains 
an excellent representation not just for the two chosen sig- 
mas of the two curves in Fig. [2l but for arbitrary values of 
a across the range spanned by these two a values. This is 
the range within which we will find the actual MDI noise 
values to lie. 

We note in Eq. (jH) that a falls between the values 2 and 
1 for the Gaussian and Lorentzian cases, as expected. 

3.2. Scaling law for the average unsigned flux density 

The unsigned vertical flux density is a resolution dependent 
quantity because the pattern of solar magnetic fields con- 
tains mixed magnetic polarities across the whole range of 
spatial scales (resolved and unresolved). Due to cancellation 
of the positive and negative contributions of the opposite- 
polartiy magnetic fluxes within each resolution element, the 
net flux or apparent flux density (as averaged over the res- 
olution element) is reduced. As we increase the spatial res- 
olution more of the mixed polarity elements get resolved, 
which reduces the cancellation effects. This causes the av- 
erage of the unsigned vertical flux density to increase. 

The scaling law that governs the variation of the av- 
erage unsigned vertical flux densit y is called the cancella- 
tion fu nction. It was introduced bv lPietarila Graham et al.l 
(|2009^ for analysis of a Hinode SOT/SP data set for the 
quiet Sun disk center , which had beeii converted to ver- 
tical flux densities by iLites et all (|2008l ) . The cancellation 
function 

Save - rf"" (5) 

describes how the average unsigned flux density Save for a 
given region on the Sun depends on the chosen resolution 
scale d. Through numer i cal sm oothing of the Hinode data 
iPietarila Graham et al.l (j2009t ) determined the value of the 
cancellation exponent k to be 0.26. The determination is 
however sensitive to the influence of instrumental noise and 
the way in which the polarimetric data are converted to flux 
densities. An independent determination of the cancellation 
exponent from quiet-sun Hinode data led to k = 0.13, a 
value half as large (Stenflo 2011,). If, as it seems, the pattern 
of solar magnetic fields behaves like a fractal, we expect 
K to be independent of scale size d over several orders of 



magnitude, including scales larger than the MDI resolution 
and unresolved small scales beyond the Hinode resolution. 

Random instrumental noise, which is spatially uncorre- 
lated from pixel to pixel, scales with where N is the 

number of pixels inside the smoothing window (simulated 
resolution element), implying a k of unity, representing a 
much steeper scaling law than that of the magnetic-field 
pattern. Noise has the effect of spuriously steepening the 
cancellation function. 

The 96 minute cadence SOHO/MDI data set of fuU disk 
magnetograms that we have used contains a mixture of 
magnetograms based on two distinctly different integration 
times, 1 min and 5 min. A keyword in the fits file header 
identifies which integration time was used for the given 
magnetogram. Since the noise characteristics are distinctly 
different for the two choices of integration time, we have di- 
vided the entire data set into two subsets, each representing 
the respective integration time. 

In their recent an a lysis of the noise in the MDI mag- 
netograms, iLiu et all (|2012l) estimated the disk-averaged 
value for the standard deviation a of the noise distribution 
to be 16.2 and 26.4 G for the magnetograms with 5-min 
and 1-min integrations, respectively. They considered these 
values to only be upper limits, since they had been ob- 
tained by fitting Gaussians to the PDF cores in quiet solar 
regions, assuming that the solar contribution to the core 
width can be neglected. If this assumption is incorrect, the 
actual noise is smaller. 

The values of lLiu et al.l ()2012f ) are for the disk-averaged 
line-of-sight component B\\ , while our analysis is based on 
the vertical flux density By = within radius vector 

r/vQ < 0.9. Since the average of l//i within this radius 
vector is 1.39, we should multiply by this factor to trans- 
late the Liu et al. values to a scale that is comparable to 
ours. Their values then become 22.5 and 36.7 G. Note how- 
ever that this crude translation only represents an estimate, 
because the noise is not constant over the disk but has a 
large-scale pattern. Nevertheless it serves as a reference to 
check the consistency of the analysis. As we will see in the 
next subsection, our use of the scaling laws to determine 
the noise a leads to the values 18.8 and 33.2 G, which are 
10-20 % smaller than the corresponding values of lLiu et"aLl 
()2012 '). As their results represent upper limits, our results 
are fully consistent with theirs. 

In Fig. [3] we illustrate for two representative magne- 
tograms how the two scaling laws, with k = I for the noise 
and K — 0.13 for the magnetic fields, combine to give a scal- 
ing behavior that can be fitted to the observational data. 
The average unsigned vertical flux density over r/rQ < 0.9 
is plotted as a function of the relative scale d, in units of 
the MDI pixel size, d represents the side of the square- 
shaped window used for numerical spatial smoothing. The 
solid curves are the actually observed values i?app of the 
average unsigned fiux density, the dotted curves the in- 
ferred intrinsic average unsigned flux density i?avo, while 
the dashed curves represent the modeled combination of 
Save with a noise a (for d = 1) of 18.8 G (upper panel) and 
33.2 G (lower panel). We see that the agreement between 
the model (dashed) and observations (solid) is nearly per- 
fect in both cases. It is through fitting of the model to the 
observations that the noise a is determined, as will be de- 
scribed next. This determination also leads to a verification 
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Fig. 3. Average unsigned vertical flux density vs. relative 
size (in pixel units) of the smoothing window, illustrated for 
two representative magnetograms, one obtained with 5-min 
integration (upper panel), the other with l-min integra- 
tion (lower panel). The solid curves represent the observed, 
noise- affected values i?app, the dotted curves the noise- free 
values Save, while the dashed curves are obtained through 
combination of i?avo with the noise a according to the model 
described by Eqs. (O and (g]). 
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Fig. 4. Standard deviation a of the Gaussian noise distribu- 
tion (vertical axis) vs. noise-free average unsigned vertical 
flux density, determined from least-squares fits to the ran- 
domly selected 320 MDI magnetograms, assuming a can- 
cellation exponent k = 0.13. The filled circles have been 
obtained from magnetograms recorded with 5-min integra- 
tion, the open circles from magnetograms with l-min inte- 
gration. The solid horizontal lines at 18.8 and 33.2 G rep- 
resent the average values of a for the two populations. For 
comparison the corresponding up per limits to the noise lev- 
els as derived from the analysis of lLiu et all (|2012[ ) are given 
as the dotted lines. 



3.3. Determination of the noise levels 

For the model fitting to determine the noise contribution 
to the observed average unsigned flux densities i?app we 
have randomly selected 320 MDI magnetograms evenly dis- 
tributed over the 15 yr period of MDI operation. The fits 
header keywords reveal that 140 of these magnetograms 
represent recordings with 5-min integration, 180 with 1- 
min integration. The magnetograms have then been spa- 
tially smoothed with square windows having side d = 3, 7, 
11, 15, 19, and 23 pixels. For each magnetogram there are 
thus 7 versions with different degrees of smoothing. They 
are all converted to vertical flux densities By through di- 
vision by ^ for each pixel. \By\ is then averaged over the 
disk region r/r© < 0.9 to give the noise-affected average 
unsigned vertical flux density iJapp as a function of d. For 
each magnetogram we thus obtain 7 observables (which in 
Fig. 13] are represented by the filled circles and solid lines, 
for the choice of two magnetograms). 

These 7 observables serve as input for an iterative least 
squares model fitting. The fit model is based on Eqs. ([3]), 
([4]), and ([5]), with the i?app values for the 7 different d values 
being the observables to be fit. For convenience we intro- 
duce the notation Bsoiar to refer to the value of Save for the 
unsmoothed case (to distinguish it from the more general 
notation Save, which may refer to any smoothed case and 
therefore depends on d) . Then the d- independent model pa- 
rameters to be determined are Ssoiar, c, and cancellation 
exponent k. For reasons of numerical robustness we only let 
.Bsoiar and a be free model parameters, while keeping the 
value of K fixed, but repeat the fit procedure for a sequence 
of K values to demonstrate that only one particular value 
of K leads to physically acceptable results. 



Let us here mention the technical detail that Ssoiar, the 
unsmoothed version of Save, actually refers to d = 2, since 
according to the sampling theorem there must be at least 
2x2 pixels per resolution element. Therefore the resolution 
does not increase when we go from d = 2 to d = 1. In 
contrast, the sampling theorem does not apply to the noise, 
which does not contain spatial structures. In Fig.|3]we have 
plotted the value of Bsoiar at d = 1 (as the left ends of the 
dotted lines), because it relates to the values of Bapp and 
(7 at d = 1. The practical consequences of the sampling 
theorem are however insignificant, since the solar scaling is 
so small when going from d = 1 to 2. Thus 2" is only 1.09 
for K = 0.13. 

Our least squares fitting gives us for each of the 320 
magnetograms the two determined fit parameters Bsoiar 
and (T, the standard deviation of the Gaussian noise dis- 
tribution. In Fig. 2] we have for a cancellation exponent 
K = 0.13 plotted cr as a function of i?soiar, using filled cir- 
cles for the magnetograms that were obtained with 5-min 
integrations, open circles for the magnetograms based on 
l-min integrations. We see that the filled and open circles 
form two distinctly different populations, as expected. The 
two horizontal lines mark the average values of cr for each 
population, 18.8 and 33.2 G, respectively. For comparison 
we give as the dott ed lines the uppe r limits to cr, 22.5 and 
36.7 G, derived bv iLiu et al.l (|2012| ) but here scaled with 
the average value of 1/fi inside r/r© < 0.9 to make them 
refer to our vertical flux density scale. They lie 10-20% 
higher than our solid lines, which is fully consistent with 
the circumstance that they represent upper limits only. 

The two magnetograms that we selected for illustra- 
tion of the scaling behavior in Fig. [3] were chosen to have 
.Bsoiar values in the middle of the range spanned by the 320 
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Fig. 5. Regression line slopes with error bars for the point 
populations in diagrams of the type of Fig. HI based on 
model fitting with different fixed values for the cancella- 
tion exponent n. The filled circles and solid line represent 
magnetograms with 5-min integration, while the open cir- 
cles and dashed line represent magnetograms with 1-min 
integration. Only zero slope is physically acceptable, which 
occurs when k = 0.13 (marked by the vertical dotted line). 
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Fig. 6. Time series of the noise-corrected disk average (over 
the disk region t/tq < 0.9) of the unsigned vertical flux 
density, smoothed with a 1-month wide time window (solid 
line) . Due to an MDI data gap the solid line is absent in the 
interval 1998.4 - 1999.0. The dashed curve uses a second- 
order polynomial in the sunspot number (also smoothed 
with a 1-month time window) to fit the solid curve. This 
polynomial fit gives a basal flux density (in the absence of 
sunspots) of 2.7 G. 



magnetograms. Thus the i?soiar value for the upper panel is 
15.6 G, for the lower panel 15.2 G. The corresponding points 
in Fig.|4]he almost exactly on the horizontal lines that rep- 
resent the average a for the respective point population. 

It is gratifying that there is no dependence of the fit- 
ted value of the noise a on the magnetic fiux level -Bsoiar 
of a given magnetogram. Any such dependence would be 
unphysical if the noise is of instrumental origin and is fully 
stochastic and additive, like photon noise, and since the 
brightness-flux correlation is very small when recorded with 
4 arcsec resolution (except in sunspots, which are expected 
to only have a minor influence on the global fiux aver- 
ages). However, this fiux invariance of a turns out to be 
unique for the chosen value of cancellation exponent k. For 
other values of k the point distribution in i?soiar - c space 
is tilted. We have repeated the fitting procedure for a se- 
quence of fixed K values, determined the slopes of the re- 
gression line fits to the point populations, and plotted these 
slopes with their error bars as a function of k in Fig.jS] The 
filled circles and solid line represent the magnetograms with 
5-min exposure, the open circles and dashed line the mag- 
netograms with 1-min exposure. We see that both sets of 
magnetograms give zero slope only for k = 0.13, while for 
all other choices of k an unphysical non-zero slope is found. 
It is this result that unambiguously leads us to the unique 
value K = 0.13. 

The same value of k was found from an analysis of the 
scaling behavior of Hinodc magnetograms for the quiet- 
sun disk center (Stenflo 2011). This consistency between 
the widely different Hinode and MDI data provides sup- 
port for the validity of our noise model, while indicating 
that K is indeed invariant over a large scale range, as ex- 
pected from fractal-like behavior of the magnetic field pat- 
tern. The circumstance that the same noise model with the 
same k, when applied to the quite different HMI data set 
for the quiet-sun disk center, also gives a noise cr that does 



not depend on the magnitude of the unsigned ffux density 
(cf. Fig. [T]), further supports the validity of our approach. 

3.4. Noise-corrected correlation with the sunspot number 

Having determined the two values of the noise a for the 
5-min and 1-min exposure magnetograms, we can now ap- 
ply Eqs. ([21) and (H)) to convert the noise-affected average 
unsigned vertical flux density Bapp to its noise-free counter- 
part -Bavo- Similar to the right panel of Fig. [T] we plot as the 
solid curve in Fig.|5]i3avo for the set of 5-min exposure mag- 
netograms, time-smoothed with a 1-month window. Like in 
Fig. [T] we model this curve in terms of a second-order poly- 
nomial in the sunspot number (also time-smoothed with 
a 1-month window) , deflned by Eq. ([T} . The resulting values 
for the fit parameters in units of G are: bo = 2.7, bi = 0.25, 
and 62 = —5.0 x 10"*. Except for the value of 60 (i?avc in 
the absence of sunspots) the fit parameters (&i and 62) are 
nearly identical to those of Fig. [1] For 69, however, there is 
a dramatic change, from 19.0 G in the noise-affected case to 
2.7 G in the noise-free case. The noise-corrected value of bo 
may be interpreted as a time-invariant, basal ffux density 
of the Sun. 

The corresponding time series of the noise-corrected 
-Bavc based on the use of the 1-min exposure magnetograms 
is nearly identical although somewhat noisier and is there- 
fore not illustrated separately here. The polynomial fit pa- 
rameters are slightly different: bo = 1.0, bi — 0.28, and 
62 = —6.5 X lO"**. This difference gives an indication of 
the uncertainty in the determination of the basal flux level. 
As the noise a for the 1-min magnetograms is almost twice 
as large (33.2 G) as for the 5-min magnetograms (18.8 G), 
we give considerably more weight to the analysis based on 
these lower-noise magnetograms. As we will see below, our 
2.7 G level found with the 5-min MDI magnetograms is fully 
consistent with the 3.0 G level that we find from analy- 
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sis of SDO/HMI magnetograms. The consistency between 
the MDI and HMI results supports the conclusion that the 
basal flux density of about 3 G is a real physical property 
of the Sun. 



HMI noise at. disk center 
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4. Basal flux density from HMI data 

The Helioseismic and Magnetic Im ager (HMI) on the 
Solar Dynamics O bservatory (SDO) (jScherrer et al.l[2012t 
ISchou et al.ll2012l) is the next-generation successor of the 
MDI instrument on SOHO and has been operating since 
April 2010. HMI delivers line-of-sight magnetograms with 
a cadence of 45 s and a spatial resolution that is four times 
higher than that of MDI, with a noise level that is better 
than that of the 5-min MDI magnetograms by more than a 
factor of 2. This substantial improvement in quality makes 
the HMI data particularly suited for a determination of the 
basal flux density level. 

The advantage of the MDI data set over that of HMI 
is that the MDI time series of 15 yr covers more than a 
full sunspot cycle, which allows an exploration of the cor- 
relation between disk-averaged unsigned flux density and 
sunspot number, as we have done in Figs. [T] and [51 The 
second-order polynomial fit to the sunspot number allowed 
us to identify the constant coefficient in the polynomial fit 
as the basal value of the disk-averaged vertical flux den- 
sity, since it represents the flux density in the absence of 
sunspots. This type of approach is not possible with the 
HMI data, since HMI so far only covers a small fraction of 
a solar cycle. 

The HMI data set on the other hand has the advantage 
of covering in great detail a cycle phase when the Sun has 
been unusually quiet. During this phase the average un- 
signed fiux densities should come closest to the minimum 
level that is given by the basal flux density. 

Like with the MDI data, the apparent basal flux density 
of the HMI data is dominated by the noise contribution. A 
determination of the basal flux density therefore has to be 
preceded by an accurate modeling and determination of the 
HMI noise. This is done in the following subsection. 

4.1. HMI noise at disk center 

For our HMI analysis we have selected 676 magnetograms, 
one for each day, starting April 30, 2010. The procedure 
for the determination of the standard deviation a of the 
Gaussian noise distribution with least squares model flt- 
ting based on the different scaling behavior of the noise 
and the magnetic fields is identical to the procedure that 
we used for MDI and which was described in great detail in 
Sect. 13.31 For each HMI magnetogram we thus generate 7 
different versions with different degrees of spatial smooth- 
ing. However, since we are not doing correlations with the 
sunspot number, there is no reason for averaging the un- 
signed flux density over most of the solar disk. Instead we 
only do the averaging over the central portion of the disk, 
within r/rQ < 0.1. Inside this disk-center region there is 
no signiflcant difference between the vertical and the line- 
of-sight directions, so no geometrical projection is needed 
to convert the line-of-sight component to a vertical compo- 
nent. With our previous notations, -Bapp is the apparent, 
noise-affected average unsigned vertical flux density within 
r/rQ < 0.1, while i?avc is its noise-free counterpart. Like 
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Fig. 7. Standard deviation a of the Gaussian noise distribu- 
tion (vertical axis) vs. noise-free average unsigned vertical 
flux density, determined from least-squares fits to the disk- 
center regions of a set of HMI magnetograms, based on the 
same fitting procedure as used for Fig. H) The average a 
(solid horizontal line) is found to be 8.0 G. For comparison 
we give as the dotted horizontal lines the upper and lower 
limits to t he H MI disk center noise cr, as determined by 
iLiu et all (l20Tl . 



for MDI, Eqs. (jS]), (jj]), and ([5]), with cancellation exponent 
K — 0.13, are applicable also for HMI, since these equations 
are not resolution dependent over the spatial scales that we 
are considering. 

The results of the least squares parameter fitting are 
shown in Fig. [71 where the noise a is plotted vs. i?soiar, the 
value of the noise-corrected i?avc for the unsmoothed ver- 
sion of the magnetogram. The average value of cr, marked 
by the horizontal solid line, is found to be 8.0 G. For com- 
parison we draw as the horizontal dotted lines the upper 
and lower limit s to the HMI noi se as derived from the care- 
ful analysis bv iLiu et all (j2012D . The upper limit of 8.5 G 
was obtained by Gaussian fitting to the core region of the 
observed PDF in quiet regions, while the lower limit of 
7.2 G represents the ideal, theoretically expected level based 
on Monte Carlo simulations of the photon noise in HMI. 
The perfect consistency of our results with these two limits 
strongly supports the validity of our noise model based on 
the use of scaling laws. 

Note that the HMI a of 8.0 G is small in comparison 
with the disk-averaged a of 18.8 and 33.2 G for the 5-min 
and 1-min MDI magnetograms. Therefore the domination 
of the noise contribution relative to the solar basal flux 
contribution will be much less severe for HMI as compared 
with the MDI data set. 

4.2. Basal flux density from the noise-corrected HMI time 
series 

With the determined value of cr = 8.0 G we can now, exactly 
like we did in Sect. 13.4] for the MDI data, use Eqs. ([3]) and 
(m to convert the noise-affected average unsigned vertical 
flux density i?app to its noise-free counterpart -Bsoiar (which 
represents Save for the unsmoothed magnetograms). Since 
we have analysed one HMI magnetogram per day, we get 
a time series, which we have plotted in Fig. [Sjfor the most 
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Fig. 8. Time series (2010.4 - 2011.6) of the noise-affected 
average unsigned vertical flux density (dashed curve) and 
its counterpart after noise removal (solid curve) , represent- 
ing the disk center region {r/rQ < 0.1). The solid curve 
never dips below the 3.0 G level, which can be interpreted 
as the basal flux density at the HMI resolution of 1 arcsec. 



quiet period of solar activity, from 2010.4 to 2011.6. The 
noise-corrected average flux density i?soiar is plotted as the 
solid curve, while the time series i?app before noise removal 
is shown as the dashed curve for comparison. 

While the value of the noise-free Bsoiar can be seen to 
vary significantly from day to day, it never drops below the 
3.0 G level marked by the dotted line, although the illus- 
trated period represents an unusually deep and extended 
quiet phase of the solar activity cycle. The 3.0 G level may 
therefore be interpreted as the basal unsigned vertical flux 
density of the Sun, which is always present, regardless of 
the level of solar activity. 

It is difficult to give a good estimate for the error bar 
in the 3.0 G HMI basal flux density level, since the error is 
most likely dominated by systematic effects (like the degree 
of validity of our noise model) rather than by formal fitting 
errors. In the case of the MDI analysis, the difference in the 
bo values of 2.7 G for the 5-min integration magnetograms 
and 1.0 G for the 1-min integration magnetograms gives an 
indication of the degree of significance of the determination, 
although we give more weight to the 5-min magnetograms 
because they are less noisy to start with. Similarly, we give 
still more weight to the HMI analysis, due to the higher 
quality and much lower noise of the HMI data set, and since 
we only use disk center data with no reference to sunspots. 

Still one may wonder if the HMI 3.0 G level could be 
much different, for instance close to zero, due to errors in 
our approach. This question will be clarified in the next sec- 
tion by comparing the appearance of a magnetogram that 
represents the basal flux level with one that has more than 
twice as much average flux density. If the basal flux density 
were much less than 3G, then the "basal magnetogram" 
would look much more empty, with much less visibility of 
the magnetic network, and the flux density histograms for 
the two magnetograms would differ much more than they 
do. From such considerations we can exclude a basal flux 
density level below about 2 G. 

The basal flux level is dependent on spatial resolution 
(but not on integration time, since the temporal evolution 



:iQll-06-S4TS3:59 




Fig. 9. Comparison of the disk-center portions (|a;|, \y\ < 
O.lro) of two HMI magnetograms, recorded at times 
2011.2086 (upper panel) and 2011.4796 (lower panel), when 
the average of the unsigned vertical flux density was 3.06 
and 6.65 G, respectively. The upper panel thus represents 
the case when the basal flux density level is approximately 
reached. The grey-scale cuts are at ±50 G in both panels. 



of the magnetograms within the integration times that we 
are dealing with is insigniflcant). The given values refer to 
the particular spatial resolution used (1 arcsec in the case 
of HMI) and scales with the resolution according to the law 
given by Eq. ([5]) with cancellation exponent k — 0.13. The 
HMI 3.0 G level therefore translates to 3.0 x 4-o i3 ^ 2.5 G 
for the 4 arcsec MDI resolution. This scaling is in nearly 
perfect agreement with the value of 2.7 G that we deter- 
mined from the analysis of the 5-min MDI magnetograms 
through the correlation between the unsigned vertical flux 
density and the sunspot number over the solar cycle. In the 
next section we will discuss the physical meaning of this 
result. 



5. Basal flux and the local dynamo 

While the average unsigned flux density is a convenient 
parameter to characterize the statistical properties of the 
magnetic pattern in terms of a single number, we need to 
examine the full pattern to understand what this number 
really means. For this purpose we illustrate in Fig. ^ the 
disk-center regions of two HMI magnetograms, both repre- 
senting the quiet Sun, but with different values of -Bsoiar, 
the noise-corrected average unsigned vertical flux density. 
For the upper panel, -Bsoiar — 3.06 G, nearly identical with 
our estimate of the basal flux density level, while for the 
lower panel it is 6.65 G, significantly larger although still 
quite small. The plotted regions extend in x (E-W direc- 
tion) and y (S-N direction) over zLQAtq with respect to disk 
center. The grey scale cuts in both plots are — 50G (dark) 
and -1-50 G (white). The upper and lower magnetograms are 
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Fig. 10. Amplitude-normalized PDFs for the two magne- 
tograms in Fig. [S] (solid line for the upper magnetogram, 
dashed line for the lower) and for the Gaussian noise dis- 
tribution with a = 8.0 G (dotted line). For comparison we 
have also plotted as the dash-dotted line the analytical rep- 
resent ation of the P DF from the Hinode quiet-sun analy- 
sis of IStenfiol ()2010[ ). after having convolved it with the 
a = 8.0 G Gaussian to make it comparable to the HMI 
PDFs. 

from times 2011.2086 and 2011.4796, respectively, and can 
be identified in the time series of Fig. |8] as a local dip and 
a local peak at these times. 

While the two magnetograms are qualitatively similar, 
in both cases characterized by a magnetic network on the 
scale of the supergranulation, they are strikingly different 
quantitatively, as expected from the difference of more than 
a factor of two of their i?soiar values. Note that before 
noise correction the apparent average unsigned flux den- 
sities were 8.1 and 10.8 G and therefore did not differ that 
much from each other (since these values are dominated by 
the noise contributions). 

5.1. Intermittent basal magnetic flux 

To make the comparison between the two magnetograms 
in Fig. ini more explicit and quantitative, we plot in Fig. [TUl 
the respective histograms of the vertical flux density, as the 
solid line for the upper magnetogram, as the dashed line 
for the lower one. These histograms (which we refer to as 
empirical probability density functions PDF, without im- 
plying any assumption concerning the underlying physical 
process), contain the noise contribution, since the noise- 
removal procedure that we applied to the average unsigned 
flux densities can only be used for this statistical ensem- 
ble average. For comparison we have therefore in Fig. [TU] 
plotted as the dotted curve the Gaussian noise distribution 
with the determined standard deviation cr of 8.0 G. It agrees 
almost perfectly in the PDF core region with the two solar 
PDFs, which only differ from it by their extended damping 
wings. 

The effect of noise on the observed PDF is a Gaussian 
convolution of the intrinsic, solar PDF. In principle one 
could therefore retrieve the noise-free solar PDF through 
deconvolution. The circumstance that the Gaussian noise 
profile almost perfectly matches the core regions of the solar 



PDFs means, however, that the core width of the solar PDF 
is insignificant in comparison with the noise broadening. 

Through forward modeling with Hinode quiet-sun data 
from February 27, 2007, it was found that an analytical rep- 
resentation in terms of a stretched exponential for the core 
region, supplemented by extended, quadratically declin- 
ing wings, gave an excellent representation of the observed 
PDF when co nvolved with the Gaussian noise distribution 
(|Stenfloll20ldh . With this analytical PDF one finds that the 
average unsigned fiux density is almost exclusively deter- 
mined by the PDF damping wings, although the PDF core 
region represents the vast majority of the pixels (the dom- 
inating filling factor). In fact, if one changes the strength 
of the wings by decreasing the damping parameter in the 
analytical representation, the average unsigned flux density 
may go down to values even far below 1 G, but the exact 
value of the lower limit depends on the details of the not 
so well known intrinsic shape of the narrow core region. 

It is thus the extended PDF wings that are the source of 
the determined (noise-free) values of the average unsigned 
fiux density, including the value of the basal flux density. 
The solid- line PDF in Fig.[TU]is representative of the mag- 
netic pattern at the basal flux density level. The difference 
in wing amplitude between the dashed and solid PDFs is 
fully consistent with the conclusion that the wing region 
is the source of the difference between their intrinsic -Bsoiar 
values, 6.65 and 3.06 G. In contrast, the apparent -Bapp val- 
ues of 10.8 and 8.2 G are dominated by the contribution 
from the noise-dominated Gaussian PDF core. 

For comparison with our HMI-based PDFs, we plot our 
analytical representation of the mentioned Hinode PDF, 
but only after it has been convolved with the Gaussian 
a = 8.0 G PDF. This convolution is done exclusively for the 
purpose of allowing us to quantitatively compare the rela- 
tive magnitudes of the different PDF wings on a common 
scale. This comparison reveals that the quiet region picked 
by Hinode contained significantly more magnetic fiux than 
our two HMI regions. The different PDFs in Fig. [10] illus- 
trate that there is a broad range of quiet-region fiux levels, 
and that this variable quiet-sun fiux can have nothing to 
do with a time invariant local dynamo. The qualitatively 
similar PDF shapes indicate that all of them are produced 
by the same type of physical processes, including the solid- 
line PDF that can be taken as representative of the basal 
fiux level. 

The extended PDF wings are a signature of intermit- 
tency, since the combination of high fiux density with low 
occurrence probability in the PDF wings implies significant 
spatial separation between such relatively rare, strong-field 
elements. The circumstance that the dominant contribution 
to the unsigned average fiux density always comes from the 
wings rather than from the PDF core implies that the mag- 
netic flux is highly intermittent, including the basal flux. 
The intermittency expresses itself in Fig. [HI in the form of 
discrete flux fragments that form a network-like pattern. 

It has been argued that while the PDF wings that repre- 
sent the intermittent network-type flux may be asymmetric 
between the positive and negative polarities, the PDF core 
(defined to refer to the internetwork) is always symmetric, 
and that a symmetric, invariant PDF core is a signature 
of a local dynamo ([Lit cs 2011). However, we have demon- 
strated here that when the noise contribution is removed 
from the PDF, the magnitude of the unsigned average fiux 
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density is almost entirely (beyond the 1 G level) determined 
by the intermittent PDF wings. 

The apparent existence of a time invariant basal flux 
level may seem to point towards an origin in terms of a 
local dynamo. The flux generated by the global dynamo 
however spreads to form a background field pattern across 
the whole surface of the Sun, and this background does not 
suddenly vanish when the sunspot number goes to zero. The 
decay and removal of flux from the solar photosphere has 
a certain time scale, which becomes longer when the flux 
level decreases, because flux removal has to do with encoun- 
ters between flux fragments of opposite polarities. Due to 
the extreme flux intermittency, the polarity encounters be- 
come increasingly rare as the flux density level goes down. 
Therefore the background flux does not have enough time 
to be fully removed before the next cycle injects new flux 
into the background. 

This implies that the global dynamo will always leave 
a residual, non-zero flux density that will be indistinguish- 
able from the basal flux density that we have determined. 
Therefore we can only say that our basal flux density rep- 
resents an upper limit to the small-scale flux density that 
may hypothetically be generated by a local dynamo. 



5.2. Hidden basal-type flux from the Hanle effect 

It has long been known from applications of the Hanle ef- 
fect that the solar photosphere is seething with a ubiquitous 
magnetic field that is tangled on scales much smaller than 
the telescope resolution and therefore does not sho w up in 
magnetogr ams that are based on the Zeeman effect (jStenfld 
Il982lll987[ ). The "Second Solar Spectrum" , the linearly po- 
larized spectrum that is exclusively produced by coherent 
scattering processes and which is the playground for the 
Hanle effect became accessib l e to sy stematic exploration in 
1994 (jStenflo fc Kelleilll996L Il997li through the introduc- 
tion of the ZIMPOL techno logy for high-precision imag - 
ing polarimetry ()Povellll995L l2cioil: iGandoriPer eralll2004[) . 
From the series of observing campaigns with ZIMPOL since 
then we have been able to notice how the general appear- 
ance of the Second Solar Spectrum has changed with the so- 
lar cycle, from an "emission- like" spectrum (in Stokes Q/I) 
full of intrinsically polarizing lines during solar minimum, 
to a mixed emission- and absorption-like spectrum during 
the maximum phase (the period extensiv ely documented i n 
the Atlas of the Second Solar Spectrum, lGandorfeill2000t ) . 
While the polarization of the molecular lines seemed to be 
nearly invariant with the cycle, the atomic lines exhibited 
large variations (Stenflo 2003). Another intriguing finding 
was that the turbulent field strength derived from the Hanle 
effect in molecular lines was much s maller than the value 
obtai ned from the Sri 4607 A line (jTruiillo Bueno et al.l 
[2OOI . 

In 2007 a synoptic program was started with ZIMPOL 
at IRSOL (Istituto Ricerche Solari Locarno) to explore the 
solar-cycle variations of the small-scal e, tangled magnetic 
field with the differential Hanle effect (jStenflo et al.lll998h 
in molecular C2 lines at 5141 A. Assuming a single- valued 
field strength and an isotropic angular distribution, an av- 
erage field of 7.4 ± 0.8 G was found for the period 2007- 
2009, with no evid ence for a significant temporal evolution 
(jKleint et al.ll201l[ ). This field therefore appears to have the 



characteristics of a basal flux at the scale where the tangled 
fleld resides. 

The value of 7.4 G seems to nicely relate to the 3.0 G 
basal level determined for the HMI spatial scale. Applying 
the d~'^ scaling law with k = 0.13, we flnd that 3.0 G be- 
comes 7.4 G if we go down to the 0.7km scale. This lies 
in the middle of the scale range for the hidden flux that 
was recently estimated (Stenflo. 2012ll based on an analy- 
sis of the magnetic energy spectrum of the quiet Sun over 
scales spanning 7 orders of magnitude. It was argued that 
the most logical location for the tangled hidden flux is just 
below the range of 10-100 km, where most of the kG type 
flux tubes reside. The tangled fields may exist all the way 
down to the magnetic diffusion limit around 10-100 m, but 
since their strength is expected to weaken as we go down in 
scale, their main contribution should come from the upper 
end of the range below a few km. 

In contrast to this modest value for the basal hidden flux 
derived from the differential Hanle effect in the optically 
thin molecular C2 lines, the Hanle effect in atomic lines re- 
veals a hidden field that is an order of magnitud e stronger. 
The most elaborate modeling so far (jTruiillo Bueno et all 
|2004[ ) of the Hanle depolarization observed in the Sri 
4607 A line has led to a field strength of 60 G for the hid- 
den field (assuming a single-valued PDF w it h an isotropic 
angular distribution). iTruiillo Bueno et"al] ()2004|) explain 
this difference in terms of a spatial correlation of the hid- 
den flux with the solar granulation. Their modeling shows 
that the C2 lines are formed exclusively in the interior of the 
granulation cells, while the Sr line has contributions from 
both the intergranular lanes and cell interiors. The implica- 
tion is that the hidden, tangled fleld is much stronger in the 
intergranular lanes than in the cell centers, which appears 
plausible, since the collapsed kG type flux also strongly cor- 
relates with the intergranular lanes (|Stenfloll201ll) . The flux 
of the small-scale tangled f ield ma.y ther efore be supplied 
by the decaying flux tubes ()Stenfloll2012|) . 

Th is interpretation has got support from ISnik et all 
(|2010[ ). who could record scattering polarization in molec- 
ular GN lines with sufficient spatial resolution to allow a 
weak statistical distinction between the intergranular lanes 
and the cell interiors. More Hanle depolariza tion was found 
i n the lanes, implying stronger fields there. IShapiro et al.l 
(1201 ll ) did detailed Hanle-effect radiative transfer model- 
ing of the ultraviolet CN lines and found them to give field 
strengths of similar large magnitudes as those found with 
the Sri 4607 A line. The circumstance that the molecular 
GN lines (and the Sr line) are optically thick, in contrast 
to the optically thin C2 lines, appears to be related to the 
difference in their behavior. To settle this issue we need 
high-resolution mapping of the Hanle depolarization effect 
in the Sr line, to determine if the field-strength contrast 
between the intergranular lanes and cell interiors really is 
as huge as the interpretation models indicate. 

Unfortunately it is difficult to find suitable atomic line 
pairs for applications of the differential Hanle effect, in 
contrast to many good choices of molecular line pairs. 
Interpretation of the Hanle effect in atomic lines therefore 
depends more on absolute than differential polarimetry as 
well as on the details of the radiative-transfer modeling, in- 
cluding the 3-D geometry of the atmosphere, the effects of 
the atomic collisions, and the small-scale thermodynamic 
structure. We need a differential approach to reduce this 
model dependence. 
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A controversial issue is also whether the strong (of or- 
der 60 G) hidden field revealed by the atomic lines varies 
with the solar cycle. The conspicuous change in the over- 
all appearance of the atomic lines in the Second Solar 
Spectrum with the solar cycle is e vidence for su bstantial 
cyclic variations of the hidden flu x (jStenflol [20031) . In con- 
trast, l^ujillQi^uenol^r^l] (j2004|) have assembled some ev- 
idence from various observations of the Sr i 4607 A line that 
the large fleld strength that they deduce from this line does 
not vary with the cycle. We need a synoptic program for 
the atomic lines to quantitatively settle this issue and de- 
termine how high the basal, time-invariant component of 
the hidden flux density really is. 

6. Concluding remarks 

The magnetic energy spectrum of the Sun extends over 
more than 7 orders of magnitude, from the global scales 
down to the magnetic diffusion scale around 10-100 m, 
where the magnetic field lines cease to be frozcn-in and 
decouple from the plasma (.Stenfio. ,2012,). Apart from the 
scales in the range 10-100 km, where most of the intermit- 
tent flux tubes formed by the convective collapse mecha- 
nism are expected to reside, there is no particular preferred 
length scale in this magneto-convective spectrum, which 
could form a basis for deciding what is local and what is 
global. The spectrum extends all the way down to the mag- 
netic diffusion limit, regardless of whether a local dynamo 
exists or not. 

The global dynamo is a source of small-scale magnetic 
structuring, but it also needs the small scales for its oper- 
ation. Much of the large-scale helicity is produced by the 
statistical left-right symmetry breaking in the ensemble of 
small-scale turbulent elements. The question whether a lo- 
cal dynamo is physically relevant for the Sun depends on 
whether the rate at which it can create new magnetic fiux 
from a seed field is competitive with the rate at which mag- 
netic flux is fed into the small scales by the turbulent shred- 
ding of flux that has been created by global dynamo pro- 
cesses. 

Analysis of the observed magnetic flux in MDI and HMI 
magnetograms shows that there is an apparently time in- 
variant average unsigned flux density of 3.0 G at the 1 arcsec 
spatial scale. This basal level represents an upper limit to 
the contribution of a hypothetical local dynamo. It cor- 
responds to 3.5 G when scaled to the Hinode 0.3 arcsec 
scale with the cancellation function that we have empir- 
ically found to be valid over this scale range. Since the 
average unsigned flux densities found in quiet solar regions 
with Hinode are much larger, most of the quiet-sun mag- 
netic flux that has been recorded with Hinode has nothing 
to do with a local dynamo. 

The average unsigned flux densities are determined al- 
most exclusively by the extended damping wings of the flux 
density PDFs, also in the 3.0 G case that represents the 
basal flux. The PDF wings are signatures of the intermit- 
tent nature of the field. A residual background fiux level on 
the Sun (with a small number density of intermittent fiux 
fragments) may survive over relatively long periods in the 
absence of solar activity and fresh input from the global 
dynamo. The 3.0 G level that we have found therefore does 
not necessarily imply a significant contribution from a lo- 
cal dynamo. An interesting question is what would happen 
to the basal-type background fiux during a Maunder-type 



minimum. If there is a local dynamo, it should maintain 
a certain flux density level indefinitely, even if the global 
dynamo were to be turned off. 

Analysis of the differential Hanle effect in optically 
thin molecular C2 fines indicates the existence of a time- 
invariant, basal-type fleld of strength 7.4 G, representing 
the region in the interior of granulation cells. Modeling of 
the Hanle effect observed in atomic lines (like the Sr i 4607 A 
line) and the optically thick CN lines in the UV gives field 
strengths that are larger by an order of magnitude. This ap- 
parent contradiction may be explained if the hidden field 
is much stronger in the intergranular lanes than in the cell 
interiors. It remains controversial how much of the strong 
hidden flux represents time-invariant, basal flux. To clarify 
the nature of these flelds we need to apply differential tech- 
niques and map the Hanle effect with a spatial resolution 
that resolves the solar granulation. 

If much of the strong hidden flux (with strength of order 
60 G or more) really represents a basal flux, then this large 
field strength cannot be directly connected via a fractal 
scaling law to the 3.0 G basal level that we have determined 
for the 1 arcsec spatial scale. While the 3 G value places a 
very tight upper limit to the possible contribution of a local 
dynamo to the spatially resolved fields, the 60 G value, if 
basal, offers plenty of room for a local dynamo. This leads 
us to a scenario, in which the relative contribution of the 
local dynamo is insignificant at all the scales that can be 
resolved by current instruments, but in which the local dy- 
namo becomes the dominant agent at the smallest scales 
in the unresolved domain, probably around and below the 
1-10 km scal e (where the hidden fiux is estimated to reside 
according to lStenfloll2012[ ). The generation of magnetic en- 
ergy by the local dynamo at these small scales would then 
raise the energy spectrum and flatten it in the domain that 
leads down to the magnetic diffusion limit. A synoptic dif- 
ferential Hanle observing program for photospheric atomic 
lines is needed to determine whether this scenario applies 
to the Sun. 
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